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ABSTRACT 

We have collected a large body of NIR {H band), UV (2000 A) and Ha measurements of 
late-type galaxies. These are used, jointly with spectral evolutionary synthesis models, to study 
the initial mass function (IMF) in the mass range m > 2 M©. For spirals (Sa-Sd), Magellanic 
irregulars (Im) and blue compact dwarfs (BCD), our determination is consistent with a Salpeter 
IMF with an upper mass cutoff M„p ^ 80 M0. The history of star formation and the amount 
of total gas (per unit mass) of galaxies are found to depend primarily on their total masses (as 
traced by the H band luminosities) and only secondarily on morphological type. The present 
star formation activity of massive spirals is up to 100 times smaller than that average over their 
lifetime, while in low mass galaxies it is comparable to or higher than that at earlier epochs. 
Dwarf galaxies have presently larger gas reservoirs per unit mass than massive spirals. The 
efficiency in transforming gas into stars and the time scale for gas depletion (~ 10 Gyrs) are 
independent of the luminosity and/or of the morphological type. These evidences are consistent 
with the idea that galaxies are coeval systems, that they evolved as closed-boxes forming stars 
following a simple, universal star formation law whose characteristic time scale is small (r ~ 
1 Gyr) in massive spirals and large (r > 10 Gyr) in low mass galaxies. A similar conclusion 
was drawn by Gavazzi & Scodeggio (1996) to explain the colour-magnitude relation of late-type 
galaxies. The consequences of this interpretation on the evolution of the star formation rate and 
of the gas density per comoving volume of the Universe with look-back time are discussed. 

Subject headings: Galaxies: general - Galaxies: evolution - Galaxies: spiral - Galaxies: ISM - Stars: 
formation -Ultraviolet: galaxies 
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1. Introduction 

The hierarchical (Cole et al. 1994) and the 
monolithic collapse (Sandage 1986) scenarios of 
galaxy formation make different predictions on the 
evolution of the primeval density perturbations 
which gave birth to galaxies. While in hierarchi- 
cal models small objects form first and large ob- 
jects are formed via successive merging of smaller 
structures, in the monolithic scenario galaxies are 
formed via a unique collapse of the original density 
perturbation without any interaction with the en- 
vironment (closed box). A wide variety of physical 
and morphological properties observed in nearby 
galaxies are expected from both models: while in a 
monolithic scenario ellipticals should have formed 
most of their stars in a time short compared to the 
collapse time in order to avoid the formation of a 
disc by gas-gas dissipation, in spirals the slower 
initial star formation activity allowed a significant 
fraction of gas to form a rotating disc (Sandage 
1986). In hierarchical models the different prop- 
erties of galaxies are the result of the initial con- 
ditions of the original merger. 

The predicted luminosity functions and density 
distributions of galaxies in the Universe should in 
principle allow to disentangle observationally be- 
tween these two models (Kauffmann ct al. 1999). 
Furthermore these models make predictions on 
the time evolution of the physical properties 
of galaxies, such as their structural parameters 
(e.g. bulge-to-disk ratio), their dynamics, their 
star formation history (thus their stellar popula- 
tions), their gaseous content and metal enriche- 
ment (Kauffmann & Chariot 1998 and references 
therein). An accurate determination of the de- 
pendence of these properties on z (e.g. Cowie et 
al. 1996) will allow us to discriminate between the 
galaxy formation scenarios. 

A detailed knowledge of the phenomenology of 
galaxies at z=0, which represents the present stage 
of galaxy evolution, is of primary importance for 
constraining evolutionary models. 
Ironically the "zero point" of galaxy evolution is 
still not sufhciently well determined, due to the 
lack of extensive observations of local galaxies. 
For example the systcmatics of their current star 
formation activity, gaseous content and metal- 
licity are still poorly known. In particular the 
UV and Ha observations which are representative 



of the current star formation activity of galaxies 
(Donas et al. 1987, Kennicutt 1998), the near-IR 
data which trace the old stellar population, and 
spectro-photomctric measurements (Zaritsky et 
al. 1994) which are necessary to derive the prop- 
erties of the underlying stellar population and the 
physical conditions of the interstellar medium, are 
still rare. 

With the ambitious aim of constructing a repre- 
sentative description of the physical properties of 
nearby galaxies, we undertook a multi-frequency 
survey of ~ 3500 optically selected galaxies span- 
ning the broadest possible range in morphological 
type (Ellipticals, Spirals, dE, Im and BCD), lumi- 
nosity (-22 < Mb < -13), and environmental con- 
ditions (clusters - isolated). Our own Ha, ncar-IR, 
millimetric and centimetric observations, together 
with data from the literature, were collected in a 
multifrcquency database (see Gavazzi, Boselli & 
Pierini, 1996). 

This database has been used to interpret the 

colour-magnitude relation (Gavazzi & Scodeggio 
1996) and the star formation activity of late-type 
galaxies (Gavazzi et al. 1998). The extension of 
Qwx near-IR survey to faint spirals, Im and and 
blue compact dwarfs (BCD), which is the subject 
of the present series of papers, combined with re- 
cent Ha observations of a set of late- type galaxies, 
and with the large body of available HI and CO ra- 
dio data makes it possible to study, in the present 
work, the phenomenology of the young stellar pop- 
ulation and of the gas content of a large sample of 
late type galaxies from Sa to Im-BCD. 

Previous works devoted at analysing the phe- 
nomenology of nearby galaxies have shown that 
the fraction of available gas and the activity of 
star formation (normalized to the visible stellar 
luminosity) increases along the Hubble sequence 
(Roberts & Haynes 1994; Kennicutt 1998). How- 
ever none of these works was focused on the role of 
the total mass in governing the evolution of galax- 
ies, an issue that we address in the present work. 
We do so using the observation that the near-IR H 
luminosity, which is dominated by the emission of 
the old main sequence and red giant stars (Bruzual 
& Chariot 1993), is a direct tracer of the dynam- 
ical mass of late-type systems (within the optical 
radius; Gavazzi et al. 1996a). Since the near- 
IR mass to light ratio was found not to vary with 
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mass for galaxies spanning the whole range in mor- 
phological type from Sa to Im and BCD, contrary 
to what was found in the optical, by normaliz- 
ing other observed quantities (e.g. the gas content 
etc.) to the near-IR luminosity we better remove 
than using optical luminosities the well known ob- 
served luminosity-luminosity or luminosity-mass 
scaling relations. Furthermore the decomposition 
of the near-IR light profiles are used to discrimi- 
nate between disc dominated and bulge dominated 
galaxies. This new approach is made possible by 
the availability of a large body of near-IR observa- 
tions of nearby galaxies described in papers I, II, 
III and IV of this series. Furthermore the available 
HI and CO data allow a precise determination of 
the total gas content for a sample of galaxies with 
a broad coverage in morphological type and lumi- 
nosity. 

The sample used and the data coverage is de- 
scribed in Section 2. Early- type (type < SOa) 
galaxies are excluded from the present analysis 
since their UV emission, being dominated by the 
extreme horizontal branch stars and their progen- 
itors (Dorman et al. 1995), does not give infor- 
mation on the young stellar population. Moreover 
we restrict the present analysis to galaxies which 
do not show signs of interaction with their sur- 
roundings. For this we use the Hi-deficiency pa- 
rameter to discriminate between "normal" galax- 
ies and galaxies suffering for gas depletion due to 
ram pressure and we include in the present study 
only isolated objects or cluster galaxies with a Hl- 
deficiency (defined as the ratio of the HI mass to 
the average HI mass of isolated objects of similar 
morphological type and linear size (Haynes & Gio- 
vanelh 1984)) < 0.3. 

The current star formation rate (SFR) and the ini- 
tial mass function (IMF) of late-type galaxies are 
studied in Section 3. The dependence of the total 
(atomic plus molecular) gas content as a function 
of morphological type and luminosity is analysed 
in Section 4. A model for the star formation his- 
tory of late-type galaxies is presented in Section 5. 
The dependence of the SFR and of the gas density 
of the Universe on look-back time is briefly dis- 
cussed in Section 6. An appendix describes how 
the Ha and UV (2000 A) data are corrected for 
[Nil] contamination and internal extinction. 



2. The sample 

Galaxies analysed in this work are taken from 
the Zwicky catalogue (CGCG, Zwicky et al. 1961- 
1968)(mpg < 15.7). They are either late-type 
(typoSOa) members of 3 nearby {vel < 8000 
km s~^) clusters (Cancer, A1367, Coma), or 
located in the relatively low-density regions of 
the Coma-A1367 supercluster (ll'^SO™ <RA< 
13^30™; 18°<dec<32°) as defined in Gavazzi et 
al. (1999a). To extend the present study to lower 
luminosities, we include in the sample the late- 
type Virgo cluster galaxies brighter than rupg 
< 14.0 listed in the Virgo Cluster Catalogue as 
cluster members (VCC, Binggeli et al. 1985). 
Furthermore VCC galaxies with 14.0 < rUpg < 
16.0 included in the "ISO" subsample described 
in Bosclli et al. (1997a) and CGCG galaxies in the 
region 12^^ <RA< 13^; 0''<dec<18° but outside 
the VCC, are considered. To avoid systematic 
environmental effects we consider the subsample 
of late-type galaxies whose Hl-deficiency is <0.3, 
typic of unperturbed, isolated galaxies. The final 
combined sample comprises 233, mainly "normal" 
galaxies (a few starburst or active galaxies might 
however be included). 

The accuracy of the morphological classification 
is excellent for the Virgo galaxies (Binggeli et al. 
1985; 1993). Because of the higher distance, the 
morphology of galaxies belonging to the other sur- 
veyed regions suffers from an uncertainty of about 
1.5 Hubble type bins. 

We assume a distance of 17 Mpc for the mem- 
bers (and possible members) of Virgo cluster A, 
22 Mpc for Virgo cluster B, 32 Mpc for objects in 
the M and W clouds (see Gavazzi et al. 1999b). 
Members of the Cancer, Coma and A1367 clus- 
ters are assumed at distances of 62.6, 86.6 and 92 
Mpc respectively. Isolated galaxies in the Coma 
supercluster are assumed at their redshift distance 
adopting Ho = 75 km s""*^ Mpc"-^. 
For the 233 optically selected galaxies, comple- 
mentary data are available in other bands as fol- 
lows: 100% have HI (1420 MHz), 99% H band 
(1.65 /xm) data available. A much sparser cover- 
age exists in the UV (2000 A) (29 %), CO (115 
GHz) (38%) and Ha (6563 A) (65%), as shown in 
Table 1 and 2. 

The morphological distribution of galaxies is 
given in Table 3. 
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2.1. Data analysis 

Hq:+[NII] fluxes obtained from imaging, aper- 
ture photometry or integrated spectra are taken 
from Kennicutt & Kent (1983), Kennicutt et al. 
(1984), Gavazziet al. (1991) Gavazzi ct al. (1998), 
Young et al. (1996) Almoznino & Brosch (1998), 
Moss et al. (1998), Heller et al. (1999) and refer- 
ences therein. Ha fluxes from Kennicutt & Kent 
(1983), Kennicutt et al. (1984) and Gavazzi et 
al. (1991) have been multiplied by 1.16, as sug- 
gested by Kennicutt et al. (1994), in order to 
account for the continuum flux overestimate due 
to inclusion of the telluric absorption band near 
6900 A in the comparison filter. Additional ob- 
servations of 66 galaxies have been recently ob- 
tained by us during several runs at the Observa- 
toire de Haute Provence (France), at San Pedro 
Martir (Mexico) and at Calar Alto (Boselli et al. 
in preparation; Gavazzi et al. in preparation). 
The estimated error on the Ha+[Nn] flux is ^ 
15%. 

The UV data (at 2000 A) are taken from the 
FOCA (MiUiard et al. 1991) and FAUST (Lamp- 
ton et al. 1990) experiments; UV magnitudes are 
from Deharveng et al. (1994) and Donas et al. 
(1987; 1990; 1995; in preparation). The estimated 
error on the UV magnitude is 0.3 mag in general, 
but it ranges from 0.2 mag for bright galaxies to 
0.5 mag for weak sources observed in frames with 
larger than average calibration uncertainties. 

HI data are taken from Scodeggio & Gavazzi 
(1993) and Hoffman et al. (1996, and references 
therein). HI fluxes are transformed into neutral 
hydrogen masses with an uncertainty of ~ 10%. 

CO data, used to estimate the molecular hydro- 
gen content, are from Boselli et al. (1997b), Boselli 
et al. (1995) and references therein. The average 
error on CO fluxes is ~ 20%; the error on the H2 
content, however, is significantly larger (and diffi- 
cult to quantify) due to the poorly known CO to 
H2 conversion factor (see Boselli et al. 1997b). 

NIR data, mostly from Nicmos3 observations, 
are taken from this series of papers: Gavazzi et al. 
(1996b,c: Paper 1 and 11), Boselh et al. (1997a), 
Gavazzi et al. (2000a, Paper HI), Boselli et al. 
(2000, Paper IV). From these data we derive to- 
tal (extrapolated to infinity) magnitudes Ht de- 
termined as described in Gavazzi et al. (2000b: 
Paper V) with typical uncertainties of ~ 10 %. 



These are converted into total luminosities using: 
IoqLh = 11.36 — OAHt + 2logD (in solar units), 
where D is the distance to the source (in Mpc). 
For a few objects we derive the H luminosity from 
K band measurements assuming an average H—K 
colour of 0.25 mag (independent of type; see Paper 
HI). 

A minority of the objects in our sample have an 
H band magnitude obtained from aperture pho- 
tometry, thus with no asymptotic extrapolation. 
For these we use the magnitude i725 determined 
as in Gavazzi & Boselli (1996) at the optical radius 
(the radius at which the B surface brightness is 25 
mag arcsec"^) which is on average 0.1 magnitudes 
fainter than Ht (Gavazzi et al. 2000a, b). 

The total H magnitudes are corrected for in- 
ternal extinction according to Gavazzi & Boselli 
(1996). No correction has been applied to galax- 
ies of type > Scd. The model independent near-IR 
concentration index parameter C31, defined as the 
ratio of the radii containing 75 % to 25 % of the 
total light of a galaxy, will be used throughout 
the paper to discriminate between disc dominated 
and bulge dominated galaxies. As shown in paper 
V, pure exponential discs are characterized by C31 
< 3, while C31 is > 3 in galaxies with prominent 
bulges. 

3. The SFR in galaxies 

Various indicators of the star formation activity 
of late-type galaxies have been proposed in the 
literature (see for a review Kennicutt 1998). In 
the present analysis we use the Ha and UV (2000 
A) luminosities, which are commonly accepted as 
the most direct indicators of the star formation 
rate (SFR). 

3.1. SFR from Ha and UV (2000 A) lumi- 
nosities 

The Ha luminosity gives a direct measure of 
the global photoionization rate of the interstellar 
medium due to high mass (m > 10 M0), young 
(< lO'^ years) O-B stars (Kennicutt 1983; 1990; 
Kennicutt et al. 1994). The total SFR can be 
determined by extrapolating the high-mass SFR 
to lower mass stars using an assumed initial mass 
function (IMF) ^(to): 
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^(m)= / km-°'dm (1) 

where M„p {Miow) is the upper (lower) mass cutoff 
and a the slope of the IMF. In the assumption 
that the SFR is constant on a time scale of some 
10'^ years, the SFR (in M yr'^) is given by the 
relation: 



SFRHa{M © yr-i) = Kncia, M^p, Mi„^)LHc.{eTg s"^ 

(2) 

where Knaict^ M^p, Miow) is the proportionality 
constant between the Ha luminosity Lna (in erg 
s"^) and the SFRna (in M yr"!). The value of 
KHa{oi, Mup, Miou)), which depends on the slope 
a and on the upper and lower mass cutoffs M„p 
and Miow of the IMF, can be determined from 
models of stellar population synthesis. Different 
values of Knaict, M^p) from the stellar population 
synthesis models of Chariot & Fall (1993) arc given 
in Table 4, all for a lower mass cutoff Miow = 0.1 
M0 and a solar metallicity. 

The UV emission of a galaxy at 2000 A is dom- 
inated by the emission of less recent (~ 10^ years) 
and massive (2< m < 5 M0) A stars (Lequeux 
1988). The UV emission becomes stationary if the 
SFR is constant over an interval of time as long as 
the life time of the emitting stars on the main se- 
quence, i.e. > 3 10^ years. Thus, assuming a star 
formation rate constant on time scales of some 10^ 
years, the rate of star formation SFRuv from UV 
luminosities at 2000 A {Luv) can be determined 
from the relation: 

S'Fi?[/y(M0yr-i) = Kuv{a, M^p, Miow)Luv{G^g s" 

(3) 

where Kuv{ct, Mup, Miow) is the proportionality 
constant between the UV luminosity Luv (in erg 
s-^A-i) and the SFRuv (in M0yr-i). Values 
of Kuv{a,Mup) at 2000 A for Miow = 0.1 M0 
and a solar metallicity kindly made available by S. 
Chariot are listed in Table 4. 

Ha and UV fluxes can thus bo used indepen- 
dently to estimate the rate of star formation in 
galaxies once they arc corrected for internal ex- 
tinction and for the contribution of the [Nil] emis- 
sion line (A 6548 A and 6584 A) which contami- 



nates the narrow band Ha photometry. The cor- 
rections applied to the present data are described 
in the Appendix. As discussed in the Appendix, 
the lack of far-IR data and of integrated spectra 
for all the sample galaxies prevents us to make 
accurate [Nil] contamination and extinction cor- 
rection for each single galaxy, forcing us to use 
"statistical" corrections. Given the broad distri- 
bution in the UV and Ha extinctions measured by 
Buat & Xu (1996) and given in Fig. 11 for each 
morphological class, we estimate an uncertainty 
in the extinction correction of a factor of ~ three 
for each single galaxy, but probably lower for an 
entire class of objects. This uncertainty, which 
might be responsible for part of the observed scat- 
ter in several star formation indicators, such as 
the birthrate parameter, should not be critical 
in the forecoming analysis since it is significantely 
smaller than the whole range in star formation ob- 
served in the sample galaxies, which spans more 
than 3 orders of magnitude. 

From now on we estimate SFRs (in M0yr~"'^) 
from UV and Ha luminosities using Kuv{cy-i Mup, Miow) 
and Kffaicf, Mup, Miow) from Table 4 for a=2.5, 
Mup=80 M0 and Af;o^=0.1 M0; for galaxies with 
both UV and Ha measurements (52 objects), the 
adopted SFR is an average of the results of the 
two methods. 

3.2. The IMF 

Since the Ha emission is due to massive (m > 
10 M0) 0-B stars, and the UV emission to mod- 
erate mass (2 < m < 5 M0) A stars, the two SFR 
determinations provide us with an indirect method 
for studying the IMF of spiral galaxies in the mass 
range m > 2 M0 (Buat et al. 1987; BeU & Ken- 
ijicutt 2000). This method can be applied only 
irthe SFR has been constant over the last 3 10^ 
years, corresponding to the life time of stars dom- 
inating the UV emission on the main sequence. 
It is well known that this assumption applies for 
"normal" unperturbed objects, such as the ones 
selected here, while it does not hold for interact- 
ing systems (Kennicutt et al. 1987) and starburst 
galaxies. 

For all galaxies in our sample with available Ha 
and UV data (52 objects), we compare in Fig. 1 
the ratio of the Ha to UV fluxes, corrected for 
extinction and [Nil] contamination as described in 
the appendix, with the values determined from the 
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stellar population synthesis models given in Table 
4. Different symbols are use for disc dominated 
(filled circles) and bulge dominated (empty circles) 
galaxies. The logarithm of the Ha to UV flux ratio 
is plotted versus the morphological type and the 
H luminosity. 

The average value log[fluxHa/ fluxUV{20m°A)] 
±0.25 is consistent with a IMF with slope a = 2.5 
and Mup = 80 M0, and differs significantly from 
the values obtained for an IMF with slopes 1.5 
or 3.5 (sec Table 4). Figure 1 a and b show that 
log[fluxHa/fluxUV{2000A)] is independent of 
the morphological type, mass and of the presence 
of a bulge. The value a = 2.5 is consistent with 
the Salpeter IMF (a = 2.35; Salpeter 1955). As 
shown in Fig. 1, the Ha to UV flux ratio is more 
sensitive to the slope of the IMF than to M„p. 

As discussed by Calzetti (1999), the Ha to 
UV flux ratio should increase with extinction in 

any fixed geometry. In the case of a slab model 
(absorbing dust and emitting stars well mixed in 
a disc), the observed Ha/UV {2800X) flux ratio 
overestimates the real value by a factor of 2 for 
typical Ha extinctions of ~ 1.1 mag (Calzetti 
1999). This value has however to be taken as un 
upper limit since it has been observed that, despite 
a factor of ~ 4 in the galactic extinction law be- 
tween 6563 A and 2000 A, UV and Ha extinctions 
are comparable because of a less efficient absorp- 
tion of the 2000 A photons dictated by the less 
concentrated distribution of the UV emitting stars 
relative to the ionizing stars inside the dusty HII 
regions (see Appendix) . We expect that variations 
in the extinction are responsible for an important 
fraction of the scatter in the Ha to UV flux ra- 
tio observed in Fig. 1 without introducing any 
important second order systematic effect (see Ap- 
pendix) . 

An a=2.5 is in agreement with the value ob- 
tained by Buat et al. (1987), who analysed 31 
galaxies with Ha and UV data, adopting a method 
similar to the one used in this work, coupled with 
an older version of the evolutionary synthesis mod- 
els. We do not confirm the trend they find with the 
morphological type, which might be contributed to 
by Hl-deficient cluster galaxies which we deliber- 
ately excluded from our analysis. 



4. The star formation history of galcixies 

4.1. The birthrate parameter b 

The birthrate parameter b is defined by Kenni- 
cutt et al. (1994) as the ratio of the current SFR to 
the average past SFR. This distance independent 
parameter b is given by: 

^ _ SFR _ SFR to (1 - R) 

< SFR >past Mstar 

where to and Mstar are the age and the total stellar 
mass of the disc and R is the fraction of gas that 
stars re-injected through stellar winds into the in- 
terstellar medium during their lifetime. While the 
SFR comes directly from the Ha and the UV lu- 
minosities, the remaining parameters must be es- 
timated indirectly. Since stars eject different frac- 
tions of gas into the ISM at different epochs of 
their life, the return parameter i? is a function 
of time and depends on the assumed IMF and 
birthrate history. Kennicutt et al. (1994) have 
shown that 90% of the returned gas is released in 
the first Gyr (more than half in the first 200 Myr) 
of a stellar generation for any assumed IMF. For 
this reason we can safely assume an instant recy- 
cling gas approximation in the determination of 
the birthrate parameter b, with a constant value 
i?=0.3, as determined by Kennicutt et al. (1994) 
for a Salpeter IMF. 

The total stellar mass Mgtar can be determined 
using the same procedure described in Kennicutt 
et al. (1994), but using H band luminosities in- 
stead of optical ones. Near-IR luminosities are di- 
rectly proportional to the total dynamical masses 
{Mtot/ Lh=^-&) at the B band 25 mag arcsec"^ 
isophotal radius, independent of morphological 
type (Gavazzi et al. 1996a). The same Mtot vs. 
Lh relationship is followed by BCDs, where how- 
ever the near-IR fiux might be contaminated by 
the emission of red supergiants; thus even in these 
objects the H luminosity can be properly used as 
an estimator of the total mass. 

Assuming to ~12 Gyrs, the birthrate parameter 
b comes directly from: 

^ SFR toil - R) 

LH{Mtot/LH)DM,o„t ^ ' 

where DMcont is the dark matter contribution to 



=1.43 
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the Mtot/Ln ratio at the optical radius, that wc 
assume to be DMcont=0-5, as in Kennicutt et al. 
(1994). As discussed by Rubin (1987), this value 
is independent on type and luminosity and thus it 
should not introduce any systematic effect in the 
determination of b. Sandage (1986) and Kenni- 
cutt et al. (1994) have shown that the birthrate 
parameter, when determined using B luminosities, 
increases along the Hubble sequence, but it is still 
not known how b scales with mass. 

The birthrate parameter (in logarithmic units) 
is plotted in Fig. 2 a and b versus the morpho- 
logical typo and the H luminosity. Fig. 2 shows a 
strong relationship between the birthrate parame- 
ter and these structural parameters, with late- type 
and dwarf galaxies (low mass objects) all charac- 
terized by a similar present and past star forma- 
tion rate {b ~ 1). There are a few galaxies with b 
> 1, i.e. with present star formation rate higher 
than averaged over the past. These arc low-mass, 
dwarf galaxies which might undergo episodes of 
star bursts. Massive spirals, on the contrary, have 
present SFRs significantly lower than in the past 
(6 ~ 0.1-0.01). The relationship between b and 
the H luminosity is however considerably clearer 
than with the Hubble type, whilst only a system- 
atic difference is observed between types < Sbc 
(6 < 0.25) and types > Scd {b > 0.25), with Sc 
spanning the whole range in b. 

If we limit our analysis to Virgo cluster galax- 
ies, for which the morphological classification is 
more accurate, we observe a less dispersed relation 
between b and the morphological type. A relation- 
ship between b and the H luminosity is however 
still observed inside any given morphological class, 
implying that part of the observed scatter in the 
b vs. type relationship is due to a stronger depen- 
dence of b on the H luminosity. 

On the other hand part of the scatter in the b 
vs. log{LH) relationship is due to the presence of 
a bulge; at any luminosity, galaxies with strong 
bulges (C3i>3; open circles) have lower birthrate 
parameters than pure exponential discs (C3i<3; 
filled circles). 

4.2. The gas content of galaxies 

It is well known that the fraction of atomic gas 

increases along the Hubble sequence (Roberts & 
Haynes 1994, and references therein), but it is 



still unknown how the total gas mass (atomic plus 
molecular) scales with the total mass of galaxies. 
The line emission of HI at 21 cm and of ^^00(1-0) 
at 2.6 mm can be used to estimate the total gas 
content (HI -|- H2) of the target galaxies. 

The molecular hydrogen content can be es- 
timated assuming a constant ratio between the 
CO (1-0) line emission and the H2 surface den- 
sity. In this work we follow Boselli et al. (1997b), 
but adopting the CO to H2 conversion factor 
X=1.0 10^0 mol cm-2 (K km s'^)'^ of Digel et 
al. (1996). For galaxies with no CO measurement, 
we assume that the molecular hydrogen content is 
10% of the HI, as estimated from isolated spiral 
galaxies by Boselli et al. (1997b). The total gas 
mass is increased by a 30% to take into account the 
helium contribution. The total gas content (nor- 
malized to the mass) of a galaxy depends strongly 
on the morphological type and on the H luminos- 
ity as shown in Fig. 3. Late-type, low mass galax- 
ies have a larger amount of gas (per unit mass) 
than early-type, massive discs. 

The relationship between log{Mgas/ Lh) and 
Hubble type has a slightly smaller dispersion if 
limited to the Virgo cluster galaxies which have 
more reliable morphological classifications. The 
relationship with the total mass is clearer than 
with the Hubble type. As for the birthrate param- 
eter, we observe a trend between log{Mgas/ Lu) 
and Lh inside each morphological class implying 
that part of the scatter in the gas vs. morphol- 
ogy relationship is due to the large scatter in the 
galactic mass within a given type. The residual of 
the log{Mgas/ Lh) vs. Lh relationship, however, 
is not correlated to the presence of a bulge. 

4.3. The relation between gas content and 
the history of star formation in galcix- 
ies 

A relationship between the SFR and the gas 
density is known to exist locally within the discs of 
nearby late-type galaxies (Kennicutt 1989). Here 
we show that this relation extends to global quan- 
tities, integrated over the whole galaxy. This is 
not obvious a priori, since it is well known that 
most of the HI gas reservoir is located outside the 
optical disc of spiral galaxies, where the star for- 
mation does not take place. 

Figure 4 shows the relationship between the his- 
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tory of star formation as traced by the birthrate 
parameter and the total gas reservoir (per unit 
galaxy mass). 

Galaxies with 6 ~ 1 have large amount of gas, 
while objects which formed most of their stars in 
the past {b ~ 0.1-0.01) almost exhausted their gas 
reservoirs, becoming quiescent. Figure 4 shows 
a segregation between pure discs (filled circles) 
and bulge galaxies (open circles), the former being 
star-forming and gas-rich, the latter more quies- 
cent and gas-poor. 

4.4. The SFE and the geis consumption 
timescale of galsixies 

The comparison of the current SFR to the gas 
content (atomic and molecular) gives an estimate 
of the efficiency of a given object to transform its 
gas reservoir into stars. The star formation effi- 
ciency (SFE) is defined as: 



relation (Kennicutt et al. 1994): 



SFE 



SFR 
M. 



(yr-^) 



(6) 



gas 



No relationship is observed between the SFE and 
the morphological type or the H luminosity (Fig 
5a and b). Pure exponential discs (filled circles) 
have, on average, higher SFEs than bulge dom- 
inated spirals (empty circles) of similar H lumi- 
nosity. If restricted to disc dominated galaxies, a 
trend is observed, with low mass galaxies having 
a lower SFE than giant discs. 

The SFE can vary within a given morphologi- 
cal class by up to a factor 10. The average value 
is 4 10~^° yrs~^. The SFE is a time indepen- 
dent quantity: it measures the instantaneous effi- 
ciency of transforming the gas into stars. Since the 
atomic gas has to go through the molecular phase 
to form stars inside molecular clouds, a more accu- 
rate measure of the SFE can be obtained if Mgas 
in eq. (6) is replaced by the gas instantaneously 
contributing to the process of star formation, thus 
by M{H2). Even using this modified definition of 
the SFE, we do not observe any relationship be- 
tween the SFE and other structural parameters, 
in agreement with Rownd & Young (1999). 

The SFE is proportional to the time-scale for 
gas depletion if the fraction of gas ejected by stars 
and recycled is taken into account. 

This time-scale, generally referred to as the 
"Roberts' time" (Roberts 1963) is given by the 



(1 -R) (1 - R) 



(7) 



where R is the returned gas fraction. As discussed 
previously, the R parameter changes with time, 
with the IMF and with the birthrate history. The 
determination of the future evolution of the disc 
is strongly related to the dependence of R on the 
birthrate history of a galaxy. Kennicutt et al. 
(1994) have studied how (1 — R)^^, the correction 
parameter for the determination of the Roberts' 
time, changes for different IMFs, star formation 
laws (i.e. the relationship between the gas surface 
density and the SFR, known as the Schmidt law) 
and star formation efficiencies, these three vari- 
ables being the most important parameters for the 
determination of the birthrate history of a galaxy. 
Their analysis have shown that (1 — i?)~^ is in the 
range 1.5-4 for most star forming discs charac- 
terized by a birthrate parameter 5=1-0.1, while it 
can be larger for rapidly evolving galaxies (6 < 
0.1) with low gas surface densities. Since the de- 
termination of i? as a function of the star forma- 
tion history of a galaxy is behind the scope of the 
present paper, we assume a standard correction 
factor {1 — R)~^ = 2.5, consistent with Kennicutt 
et al. (1994). 

The Roberts' time is independent of the mass 
and of the morphological type (Fig. 6). This result 
is in contrast with that obtained by Donas et al. 
(1987), who observed a weak trend between the 
Roberts' time and the morphological type, with 
late- type systems showing a longer tr. The time- 
scale for gas depletion is in the range between ~ 1 
Gyr and 25 Gyr, with an average value of 10 Gyr, 
consistent with the value found by Kennicutt et al. 
(1994). These values are upper limits to the real 
gas consumption time scales because a large part 
of the available HI gas is located outside the op- 
tical disc of galaxies and thus does not contribute 
directly to the star formation, unless gas inflow 
toward the centre of the galaxies is invoked. 

5. A model of the star formation history 

Let's summarize the evidences collected so far: 
1) photometric Ha and UV (2000 A) measure- 
ments of 52 late-type (Sa-Im-BCD) galaxies, com- 
bined with spectral evolutionary synthesis models 
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show that the IMF of galaxies is consistent with 
a Salpeter IMF of slope a = 2.35 and an upper 
mass cut off of about 80 M0, in the mass range 
m > 2 M0. Galaxies of different morphology and 
luminosity seem to have similar IMFs. 

2) the birthrate parameter h, which gives the 
present-to-past SFR ratio {SFR/ < SFR >past), 
is strongly correlated with the total mass of galax- 
ies as traced by their near-IR luminosity; the rela- 
tionship between b and the morphological type is 
weaker than with Lh- 

3) the total gas reservoir per unit mass anti- 
correlates with the total mass of galaxies; as for b, 
the relationships between Mgas/ Lh and the mor- 
phological type are more dispersed than that with 
Lh- Low mass, dwarf galaxies, which are gen- 
erally pure exponential discs, have a higher gas 
content per unit mass and present to past SFRs 
than early-type, massive spirals. 

4) the birthrate parameter b correlates with the 
total gas content of galaxies. 

5) the SFE and the time scale for gas depletion do 
not strongly correlate with properties of galaxies 
such as the luminosity and/or the morphological 
type. If limited to pure discs, however, a trend be- 
tween SFE and the total galaxy mass is observed. 
The time scale for gas depletion is ~ 10 Gyrs for 
normal, unperturbed late-type galaxies. 

Point 1 is in agreement with the universal- 
ity of the IMF in galaxies found in the high- 
mass range of normal galaxies (Scalo 1986, 1998; 
Massey 1998). Recent results based on star counts 
in OB associations in the Magellanic clouds seem 
to indicate that the IMF for massive stars is in- 
dependent of metallicity (Hill et al. 1994). No 
systematic differences have been observed in the 
IMFs of massive stars associations in the Milky 
Way, LMC and SMC (Hill et al. 1994; Massey et 
al. 1995: Hunter et al. 1996c, 1997), M31 (Hunter 
et al. 1996a) and M33 (Hunter et al. 1996b); in 
all these nearby galaxies the IMF has a slope con- 
sistent with a=2.5 for masses > 1 M0. 

Points 2 and 3 have strong implications on the 
models of galaxy formation and evolution. The 
importance of the h parameter resides on the fact 
that it gives directly an idea of the history of 
star formation of galaxies. In galaxies with a 
very small value of b (generally early-type, mas- 
sive galaxies) most of the stars have been formed 



at early epochs, and the present rate of star for- 
mation is lower than the past one. This is con- 
sistent with the idea that a rapid collapse might 
have induced a strong starburst, which efficiently 
transformed most of the gas into stars. The lack 
of gas at the present epoch makes these galaxies 
quiescent. Conversely, in objects with 6^1 (late- 
type, dwarf, low mass galaxies) the gas is presently 
transformed into stars at the same rate than in the 
past. These observational evidences are consis- 
tent with the model of galaxy formation discussed 
by Bandage (1986), who proposed that galaxies 
are coeval systems, formed from the collapse of a 
primordial gas cloud, with a collapse time-scale 
depending on angular momentum. In elliptical 
galaxies the collapse was efficient enough to trans- 
form most of the gas into stars within a few 10* 
yr. In late- type systems the initial SFR was com- 
parable with the present one (few solar masses per 
year), so that ~ 10 Gyrs after the formation of the 
primeval galaxy a large fraction of the gas is still 
available to feed new stars. This idea was rc;ccntly 
rediscussed by Gavazzi et al. (1996a), Gavazzi 
& Scodeggio (1996), and Gavazzi et al. (1998). 
They observed that the U - B, B - V ., UV - B 
and B — H colour indices and the Ha equivalent 
width depend strongly on the galaxy mass. They 
interpreted these observational evidences as an in- 
dication that the total mass contributes to the 
regulation of the process of collapse of primeval 
gas clouds. This idea is also supported by the re- 
cent numerical simulation by Noguchi (1999), in- 
dicating that the bulge-to-disc ratio seems primar- 
ily regulated by the total mass of the galaxy. A 
similar result was obtaind by Boissier & Prantzos 
(2000) and Boissier et al. (2000), whose chemo- 
spectrophotometric models of galaxy evolution re- 
produce several observed properties of disc galax- 
ies only when the gas infall (and thus the star for- 
mation activity) is regulated by the total mass of 
the galaxy. 

In the assumption that galaxies are coeval, 
and that they evolved as closed boxes, Gavazzi 
& Scodeggio (1996) and Gavazzi et al. (1998) 
were able to reproduce the relationships between 
the colour indices, the Ha equivalent widths and 
the H luminosities, using the population synthesis 
models of Bruzual & Chariot (1993) and Kenni- 
cutt et al. (1994) and adopting a star formation 
law of the type: 
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SFR{t) = SFRoC 



(-t/r) 



M /yr (8) 



with a Salpeter IMF, where t is the age of the 
galaxy, r is the time scale for collapse, and SFRg 
is the rate of star formation at the time of the 
galaxy formation (i=0). They found small values 
of T (r ~ 0.5 Gyr) for giant discs and large (r ~ 
10 Gyr) for low mass objects. 

The observational results of the present work 
are consistent with this simple picture. Let us as- 
sume that galaxies evolve as closed boxes following 
a star formation law of the type described in eq. 
(8), but assuming a more realistic age of to = 12 
Gyrs. 

Mstar{t) at a given time t is given by: 



If Mtot is the total mass of a galaxy (constant with 
time in a closed box scenario), we can assume that: 

Mtot = Mgas (t) + Mstar {t) + Mdm (t) Vi 

(12) 

where Mgas , Mstar and M]j m arc respectively the 
total gas, star and dark matter masses at a time 
t. If we assume that all the gas will be in the 
future transformed into stars, for t ^ oo Mgas{t) 
0, while Mstar{t) = for i=0, thus Mstarioo) = 
MgasiO) = Mtot - Mdm{0) = SFRoT (assuming 
that all the returned gas is used to form stars). If 
we assume that Mdm does not change with time, 
Mgas{t) can be determined from relation (12). 
By substituting SFRoT to Mtot 
becomes: 



Mdm, eq. (12) 



M, 



(i) = {1-R) I SFRoc'^-^'/'^Ut' = (1 - i?)5Fi?oT[]anye(-*/^)] 
Jo 

(9) 

^ that can be directly measured from the H band 
luminosity as described in sect. 4.1: 



Mgas{t) = SFRoTB^-''^^ 



(13) 



^^-(*)/^^(*)= (l-i,)(l-eM/.)) (14) 



Mstarit) = LH{t)DM,ont(t)[Mtot/LH{t)] (10) 

given the fact that the H luminosity traces the 
total dynamical mass of galaxies. SFRo gives the 
number of stars formed per year at the formation 
of the galaxy. Its dependence on the total mass 
implies that the process of collapse is not scale- 
free. To simplify the formalism, we define X{t) = 
DMcont{i){Mtot/ Lh)- As discussed in sec. 3.2, 
X(i)=0.5 X 4.6 at t = io=12 Gyr. 
Combining eq. (5) with eq. (8) and (9) we can 
directly estimate the birthrate parameter 6 as a 
function of t and r: 



in the assumption that R does not change with 
time. 

At the same time this simple model predicts 
that the SFE defined in eq. (6) should be: 



SFE 



(15) 



m = 



t e(-*/") 
-(1 -e(-*/^)) 



(11) 



The values of b obtained from eq. (11) for various 
values of r and assuming f=to=12 Gyrs are con- 
sistent with those obtained from the population 
synthesis models of Kennicutt et al. (1994). 



thus independent on t. This means that each 

galaxy should have had an efficiency in transform- 
ing gas into stars which depends only on the r 
parameter, but similar at any epoch. 

This analytic representation of all the observed 
variables used in this work {SFR, b, Mgas, Lh, 
SFE) can be done once a star formation law such 
as that given in eq. (8) is assumed ^. 

The analytical model can be compared to the 
observables only once SFRo is known, Lh be- 
ing a function of SFRo; the other variables (6, 



^as discussed in sect. 4.1 R can be taken as a constant for t 
> 1 Gyr 



^We liave repeated this exercize replacing the exponential 
star formation law (eq. 8) with a "delayed exponential" 
law, such as the one proposed by Sandage (1986). This law 
better accounts for cases with increasing SFR with time, as 
observed in some irregular galaxies such as IZwlS (Kunth 
& Ostlin 2000). We found results qualitatively consistent 
with the exponential case, thus we avoid to expand the 
"delayed exponential" case in full details 
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Mgas/Ln, SFE), being normalized entities, do 
not depend on SFRo- This parameter, which gives 
the rate of star formation at the time of the galaxy 
formation (t=0) , is clearly dependent on the mass 
of the galaxy {SFRo xt = M.tar + Mgas). 

Following Gavazzi & Scodcggio (1996), if we as- 
sume a simple empirical relationship between the 
H luminosity and the exponential decay time scale 
for star formation of the type: 



logLff = a X logr + c 



(16) 



we can derive from eq. (9) an (10) a semi- 
analytical relationship between SFRg and r: 



SFRo 



X(fo)T°10° 

(1 - i?)re(-i2/-) 



(17) 



Once SFRo is determined, we can compare 
the predictions of the analytical models with 
the relationships previously discussed between b, 
Mgas/Ln or SFE and the H band luminosity 
(Fig 7 a,b,d). The model prediction (dotted line) 
can be directly compared to the observations for 
the Mgas/Lu vs. b relationship (Fig. 7 c), with- 
out any assumption on the relationship between r 
and Ljj, both being normalized entities thus inde- 
pendent of SFRo. At the same time the obvious 
observed scaling relationships between SFR, Mgas 
and Lh can be compared to the models prediction 
(Fig. 7 e,f) and used, together with the normal- 
ized entities, to constrain the free parameters a 
and c. An excellent fit to the data is obtained for 
a = -2.5 and c = 12 once r is expressed in Gyrs 
consistent with Gavazzi & Scodeggio (1996; a = 
-2.5; c=11.12). 

In conclusion, the observed relationships be- 
tween optical, near-IR, UV colours, Ha E.W.s, the 
birthrate parameter &, the total gas mass and the 
H luminosity can, at least to the first order, be 
reproduced in the simple assumption that galax- 
ies arc coeval (r-^ 12 Gyrs), that they evolved as 
closed boxes following an exponentially declining 
star formation law, with a decay time scale de- 
pending on the mass of the primeval cloud, small 
(r ~ 0.5 Gyrs) for massive objects {Lh = 10^^ 
Lhq), and large (r ~ 10 Gyrs) for dwarf systems 
{Lh = W^Lhq). 

If the models here discussed are realistic, we 
conclude that the strong relationship observed be- 
tween various star formation tracers, such as the 



UV flux (Donas ct al. 1990; Dcharvcng ct al. 1994; 
Buat 1992; Boselh 1994), the Ha flux (Kennicutt 
1989; 1998; Scodeggio k Gavazzi 1993 and ref- 
erences therein) and the HI or total gas surface 
density in unresolved galaxies is not due to a lo- 
cal relationship between the gas column density 
and the process of star formation, known as the 
Schmidt law, but it is simply a consequence of the 
closed-box evolution. From eq. (8) and (13) we 
expect in fact that the integrated SFR and the to- 
tal gas content decrease with time with a similar 
decline. This conclusion is corroborated by the 
fact that in "normal", unperturbed galaxies the 
HI gas (which constitutes ~ 90% of the total gas 
reservoir; Boselli et al. 1997b) is found in a flat 
disc of diameter about twice as large as the optical 
one (Cayatte et al. 1990). This gas is not expected 
to contribute directly to the star formation since i) 
up to 75% of the HI reservoir is located outside the 
optical disc where the star formation takes place; 
ii) the HI has to go through the molecular phase 
before forming stars. Thus an infall of gas from the 
outer HI disc to the optical disc must be present. 
Spiral galaxies of all morphological types and lu- 
minosity still have enough gas to continue their 
present day star formation activity for ~ 10 Gyrs. 
This time is significantely longer than previous es- 
timates because we assume an higher contribution 
of the recycled gas (see Kennicutt et al. 1994) 
and because the Hl-deficient cluster galaxies are 
excluded from the analysis (Boselli 1998). 

6. Comparison with observations at high 
redshift 

The aim of this section is to discuss if the 
closed-box scenario discussed in the present pa- 
per, based on observations of local "adult" late- 
type galaxies, is consistent with observations at 
larger redshift. 

Several recent attempts to reconstruct the evo- 
lution with look-back time of the total star forma- 
tion activity of the Universe (Cowic ct al. 1999; 
Steidel ct al. 1999; Madau ct al. 1998 and refer- 
ences therein) unanimously concluded that it in- 
creases by about an order of magnitude from z = 
to 2 ~ 1. The dependence for z > 1 is more con- 
troversial, however the most recent determination 
by Steidel and collaborators (1999) indicates that, 
if the data are appropriately corrected for extinc- 
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tion, the star formation activity of the Universe re- 
mains roughly constant for larger look-back times, 
up to 2; ~ 4. 

Eq. (8), together with eq. (16) give the time 
evolution of the SFR of galaxies of luminosity Lh, 
for a given SFRo- 

We can predict, for a galaxy of a given H lumi- 
nosity at z=0, its SFR at any z. We can then es- 
timate the SFR per unit covolume of the Universe 
by integrating the contribution to the SFR of each 
galaxy, weighted according to a given luminos- 
ity function. Since we are dealing with late-type 
galaxies, we assume the local B band Schechter 
luminosity function of Heyl et al. (1997) for types 
Sa-Sm and transform B magnitudes into H magni- 
tudes assuming the B — H vs. B colour-magnitude 
relation: 

B-ff = 7.7-0.38x5 (18) 

as determined from the subsample of Virgo late- 
type galaxies which spans the full dynamic range 
in luminosity. 

The resulting SFR per comoving volume of the 
Universe is compared with the observed SFR as a 
function of z in Fig. 8. For consistency with previ- 
ous works we assume A=0, Ho=50 km s~^ Mpc""'^ 
and go=0.5. The "observed" values of SFR (in 
M yr-iMpc-3) ^^^^ Steidel et al. (1999) ( open 
symbols) are corrected for extinction. The values 
of the 2800 A luminosity density given by Cowie 
et al. (1999) (filled squares) are transformed into 
M yr~-'^Mpc~'^ and corrected for extinction con- 
sistently with Steidel et al. (1999). The value of 
Treyer et al. (1998) (filled triangle) and Gallego et 
al. (1996; filled circle) obtained respectively from 
the local UV at 2000 A and Ha luminosity func- 
tions, are also transformed into M0yr~^Mpc~'^ 
as described in sect. 3, adopting for simplicity an 
extinction correction of 0.3 mag. Such a small UV 
and Ha extinction derives from the assumption 
that the UV and Ha luminosity functions are dom- 
inated by low mass, low luminosity, blue galaxies 
(see Appendix). 

Figure 8 shows that, at least qualitatively, the 
model predictions give an excellent fit to the up- 
dated estimates of the SFR per comoving volume. 

Our model (dotted line) is in excellent agree- 
ment with the points at z=0 by Gallego et al. 
(1996; filled circle) and by Treyer et al.(1998). It 



accurately reproduces the relation SFR Mpc ^ a 
{1+zy-^ observed by Cowie et al. (1999) in the 
range < 2 < 1.5 (filled squares), and is con- 
sistent with the extinction-corrected estimates of 
Steidel et al. (1999) for 3 < .z < 4.5. 

It is remarkable that the predicted SFRs are 
in agreement with the observations, even though 
we consider only the contribution of the late-type 
galaxies to the evolution of the SFR of the Uni- 
verse. One should remember however that the 
adopted SFR per comoving volume, with which 
we compare our model predictions, have all been 
estimated from UV-selected samples, thus biased 
towards low-extinction galaxies. Elliptical galax- 
ies, which are expected to be formed in violent 
starbursts with strong dust extinctions, should not 
be present in UV-selected samples. Their contri- 
bution at early epochs should be however observ- 
able in the millimetric domain (Franceschini et al. 
1994). 

Using a similar procedure we can compare the 
model predictions with the observed comoving 
density of the g determined from observa- 

tions of damped Lya systems. As above, the con- 
tribution to the total gas density of the Universe 
is weighted according to the B luminosity func- 
tion of Heyl et al. (1997), and B magnitudes are 
transformed into H magnitudes using eq. (18). 
The T parameter in eq. (13) is constrained by 
the relationship between r and the H luminosity 
given in eq. (16) (with a=-2.5 and c=12). The gas 
density of the Universe at different epochs can be 
estimated by integrating the derived "gas mass" 
luminosity function. 

The models predictions are found in remark- 
able agreement (see Fig. 9) with the gas density of 
the Universe Slg (expressed in units of the present 
critical density) as determined from the statisti- 
cal analysis of damped Lya systems by Pei et al. 
(1999) and Rao & Turnshek (2000). We use the 
values of ilg given by Pei ct al. (1999) corrected 
to take into account the missing damped Lya sys- 
tems for dust obscuration in an optically selected 
sample of quasars. To be consistent with Pei et al. 
(1999) we removed the contribution of H2 (10%) to 
the values of Qg estimated from our model. Both 
models predict a gas density of the Universe at 
z=0 ~ 4 times higher than the observed value de- 
termined from the local HI luminosity function of 
Zwaan et al. (1997); this is due to the fact that the 
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models overestimate the gas content of low lumi- 
nosity galaxies (see Fig. 7f), whose contribution 
to the local luminosity function is dominant. 

The agreement between fig predicted by our 
model and the result of the statistical analysis of 
damped Lyman alpha systems is good at redshifts 
< 2, while it becomes poor at higher z. The de- 
termination oi fig at z > 1.5 strongly depends on 
the adopted correction for missing absorbers at 
high redshifts, and thus, as for the SFR, on the 
variation of the dust to gas ratio with z. The dis- 
agreement between the model prediction and Qg 
for z > 2 could thus be due to an underestimate 
of the extinction bias at high redshifts in the op- 
tically selected sample of quasars ^. 
Pei et al. (1999) interpret the decrease of fig with 
z as due to the fact that damped Lyman alpha 
systems are galaxies which grow by accreation of 
(essentially ionized) gas. This would result into a 
decrease of the co-moving neutral gas density of 
the Universe with redshift. Our closed-box model 
cannot reproduce any decrease of ilg with z since 
the gas is always assumed in the neutral phase. 

We remark that the model predictions at high 
redshift depend strongly on the parameters a and 
c of eq.(16) and on the assumed B luminosity 
function. Furthermore the absolute value of the 
SFR Mpc~^ at various z might change by a factor 
of ^ 2 when adopting different population synthe- 
sis models. 

To summarize, this work has shown the impor- 
tance of mass in parametrizing most of the phys- 
ical properties of nearby galaxies. The history of 
star formation and the amount of total gas (per 
unit mass) of galaxies are found to depend primar- 
ily on their total mass and only secondly on their 
morphological type. These strong observational 
evidences must be reproduced by the models of 
galaxy formation and evolution. 

The analysis carried out in the present paper 
can be considered as an exercize aimed at show- 
ing that the "monolithic" scenario of galaxy for- 
mation is not yet observationally falsified, both 
on local and on cosmological scales. Several scal- 
ing relations observed in nearby galaxies (i.e. the 

^We notice that the extinction correction assumed by Pei et 
al. (1999) to correct (ig for missing absorbers is lower than 
that used by Steidel et al. (1999) to estimate the SFR at 
high z. 



gas content, the activity of star formation and the 
stellar populations) can be reproduced assuming 
that galaxies are coeval and that they evolved as 
closed boxes with an universal star formation law 
whose characteristic time scale parameter r scales 
inversely with the total mass of the parent galax- 
ies. The same scenario predicts a look-back time 
dependence of the integrated star formation ac- 
tivity and gas content of the Universe consistent 
with the observations up to large redshifts. In 
order to reproduce the observed local properties 
of galaxies and the cosmological dependences dis- 
cussed above, we see no compelling need for more 
sophisticated galaxy formation models, such as the 
hierachical model proposed by Kauffmann et al. 
(1993) and by Cole et al. (1994). 

The present work seems to indicate that the 
evolution of the star formation rate and the gas 
density of the Universe is due to the gas consump- 
tion via star formation processes according to a 
simple, universal star formation law, independent 
on z, valid for all rotating systems. 

Beside any cosmological speculation, the firm 
conclusion of the present work is that the history 
of star formation of disc galaxies is primarily reg- 
ulated by their total mass and only secondly by 
their angular momentum. 
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A. Corrections to Ha+[NII] and UV fluxes 
A.l. The [Nil] contamination 

The analysis of the integrated spectra of late- type galaxies by Kennicutt (1992) has shown that, on 
average, the [NII]{6583A)/Ha emission line ratio is ~ 0.5. The sample of Kennicutt (1992), however, 
is dominated by active galaxies (see sect. A. 2). In order to estimate a [Nil]/ Ha ratio representative of 
normal galaxies, as those analysed in the present work, we plot in Fig. 10 the [Nil]/ Ha ratio for spirals 
and irregular galaxies in the sample of Kennicutt (1992) as a function of the Hubble type, excluding peculiar 
objects such as Markarian galaxies. 

The figure clearly shows a different trend between spirals and irregulars, with an average [Nil]/ Ha = 
0.42 ± 0.19 for Sa< type < Scd and [NII]/Ha = 0.25 ± 0.15 for > Sd. These values are consistent with the 
result of Kennicutt (1983) based on high resolution spectroscopy of bright HII regions in nearby galaxies: 
\NII]/Ha = 0.33 ± 0.12 for spirals and [Nil]/ Ha — 0.08 ± 0.05 for irregulars. In the present analysis we 
correct the Ha fluxes for [Nil] contamination using [Nil]/ Ha = 0.42 for galaxies of type <Scd (including 
Peculiar galaxies and unclassified spirals) and [Nil]/ Ha — 0.25 for galaxies of type > Sd. 

A. 2. The UV and Ha extinction in normal galaxies 

Ha and UV fluxes can be used to estimate the rate of star formation in galaxies provided that they are 
corrected for internal extinction. Independent works based on different techniques indicate that the average 
extinction of the Ha and UV (2000 A) fluxes of normal galaxies are similar, despite a difference by a factor 
~ 4 in the extinction coefflcient kx, as determined from the galactic extinction law (Savage & Mathis 1979; 
Bouchet et al. 1985). This observational evidence is generally explained as a result of a combined effect of 
a standard extinction curve with a different geometrical distribution of emitting stars and absorbing dust 
at 2000 A and in Ha. While the Ha emission of late-type galaxies is dominated by the photoionization of 
gas by hard UV photons produced by massive (to >10 M0) 0-B stars (Kennicutt 1998) located inside HII 
regions, the UV emission at 2000 A is dominated by A stars of intermediate mass (2 < m < 5 M0) (see sect. 
3), old enough to have migrated outside the HII regions. This has been observed for example in M33 (Buat 
et al. 1994), where ~ 20% of the total UV emission at 2000 A is diffuse, even though part of the diffuse 
emission might be scattered light. Given the complex distribution of absorbing dust and emitting stars, 
the determination of the fraction of scattered light in spiral galaxies is very uncertain; however the models 
of Witt & Gordon (2000) indicate that the scattered light in the UV at 2000 A is ~ 20% once a clumpy 
geometry is assumed. Furthermore the observations of nearby galaxies are limited to few late-type spirals, 
preventing us to quantify the contribution of the diffuse emission to the total UV emission in early-type 
spirals. It is however conceivable that in Sa the diffuse emission contribution is even higher than in Scd 
galaxies such as M33. Since dust is mostly associated with HII regions, the relative geometrical distribution 
of absorbing dust and emitting stars, which increases the absorption of Ha with respect to UV photons at 
2000 A, should compensate the ~ 4 times higher UV extinction with respect to extinction at Ha. 

Prom the analysis of the Balmer decrement of a large sample of galaxies with long-slit integrated spec- 
troscopy. Kennicutt (1992) estimates that the average extinction of the Ha flux of late-type galaxies is ^ 1 
mag. Consistently the comparison of free-free radio fluxes and Ha fluxes of normal galaxies by Kennicutt 
(1983) indicates an extinction of 1.1 mag. In order to study whether the Ha extinction is morphological-type 
dependent, we have re-examined the sample of Kennicutt (1992) limited to "normal" galaxies (thus excluding 
Markarian objects), dividing galaxies in two subsamples: Sa< type < Scd and type > Sd. High resolution 
spectra were selected as first priority. The extinction of the Ha flux strongly depends on the assumption 
on the underlying Balmer absorption, in particular in normal galaxies where the H/3 line is generally weak. 
Values found in the literature range between 2 and 5 A in H(3E.W. (Kennicutt 1992; McCall et al. 1985). 
We calculate for simplicity the Ha extinction using these two values of the underlying Balmer absorption. 
Including this correction and following the prescription of Lequeux et al. (1981) we can estimate that the 
extinction at a given A from the Balmer decrement is: 
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A{X) = 2.5C(A) (Al) 



where 



C{\) 



1 r I{Ha) 



log 



F{Ha) 



+ log{l + 



E 



) 



(A2) 



HI3E.W. 



where f{Ha) is the galactic reddening function normahzed at H/3, f{Ha)=-0.355 (Lequeux et al. 1981), 
I{Ha) and I{H(}) are the extinction corrected intensities and F{Ha) and F{H(5) are the observed intensities 
of the Ha and H/3 lines. E is the underlying Balmer absorption (in A). Assuming a ratio of I{Ha)/I{HP) 
=2.86 for a case B nebula at 10'' K (Brocklehust 1971), we obtain an average value of the Ha extinction of 
A{Ha) = 0.78 ± 0.47 (for E = 5 A) and A{Ha) = 1.60 ± 0.77 (for E = 2 A) for spirals, and A{Ha) = 0.41 
± 0.29 (for £ = 5 A) and A{Ha) = 0.86 ± 0.42 (for = 2 A) for galaxies with type > Sd (see Fig. 11). 

We did not observe any relation of A{Ha) or Ha surface brightness (extended to all the galaxies of our 
sample) with inclination, even once galaxies are divided into different morphological class bins. We stress 
however the fact that the sample of Kennicutt (1992) is biased towards strongly star forming galaxies, even 
if Markarian galaxies are excluded, and thus it could be not representative of the sample of galaxies analysed 
in this work. In fact the average Ha + [NIIjE.W. in the sample of Kennicutt (1992) is 82 A for spirals and 
105 A for types > Sd, while in the sample analysed in this work the average values are Ha + [NII]E.W. 
= 22 A and Ha + [N I I]E .W. = 46 A respectively. In conclusion we decided to assume a standard Ha 
extinction correction of 1.1 mag for Sa< type < Scd, and 0.6 mag for type > Sd. 

The lack of any relationship between the UV surface brightness and inclination observed for our sample 
galaxies prevents us to estimate an inclination-dependent extinction correction law. As shown by Buat & 
Xu (1996), the UV extinction at 2000 A is ~ 0.9 mag for galaxies with type < Scd and ~ 0.2 for galaxies 
with type > Sd. If however the ratio between the Ha (corrected for [Nil] contamination) and the UV fluxes 
for galaxies with type Sa-Scd is plotted versus their axial ratio, a residual trend is still marginally present 
(see Fig. 12). 

A linear fit to this relation gives; 



where log{b/a) is the axial ratio. If we make the reasonable assumption that this trend is due to increasing 
extinction with the inclination, this relation can be used to determine an empirical correction to the UV data 
as a function of the axial ratio. This correction, however, represents a lower limit since it has been determined 
assuming that the Ha extinction is independent of the galaxy inclination. For a random distribution of (b/a) 
we derive an average UV extinction of A{UV)i ~ 0.3 mag. In order to be consistent with Buat & Xu (1996), 
we assume: 



where the values of A{UV) f.o and ki are given in Table 5 for different morphological classes. Given the 
systematic difference in luminosity between spirals and dwarfs, the adopted corrections takes into account, 
at least to the first order, a possible dependence between extinction and luminosity. 



log{LHa/Luv) = -0.56 x log{b/a) + 1.24 



(A3) 



A{UV) = A{UV)f.o{ty) - h{ty) x log{b/a) 



(A4) 
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Fig. 1. — The logarithm of the ratio of the Ha to 
UV (2000 A) fliixes versus a) the morphological 
type, b) the H luminosity. To avoid overplotting, 
a random number between -0.4 and 0.4 has been 
added to each numerical type. Filled circles are 
for disc dominated galaxies (C31 < 3), open circles 
for bulge dominated galaxies (C31 > 3). Ha and 
UV fluxes are corrected for extinction and [NH] 
contamination. 

Fig. 2. — The relationship between the logarithm 
of the birthrate parameter h and a) the morpho- 
logical type, b) the H luminosity. Symbols as in 
Figl. 

Fig. 3. — The relationship between the logarithm 
of the total gas mass (normalized to the H lumi- 
nosity) and a) the morphological type, b) the H 
luminosity. Symbols as in Fig 1. 

Fig. 4. — The relationship between the logarithm 
of the birthrate parameter h and the total gas mass 
normalized to the H luminosity (in solar units). 
Symbols as in Fig 1. 

Fig. 5. — The relationship between the logarithm 
of the SFE and a) the morphological type, b) the 
H luminosity. Symbols as in Fig 1. 

Fig. 6. — The distribution of the time for gas de- 
pletion (Roberts' time) in 6 classes of morpholog- 
ical type. 

Fig. 7. — The comparison between the model 
prediction and the observations for the a) 6 vs. 
login , b) logMgas/Ln vs. login, c) b vs. 
logMgas/ln, d) SFE vs. login, e) SFR vs. 
login and f) Mg^s vs. login relationships. Sym- 
bols as in Fig 1. 

Fig. 8. — The relationship between the SFR per 

comoving volume of the Universe and z compared 
to the model predictions (dotted line); the ex- 
tinction corrected values are from: Lilly et al. 
(1996) [open circles], Connolly et al. (1997) [open 
squares], Madau et al. (1997) [open triangle], Stei- 
del et al. (1999) [open stars], Treyer et al. (1998) 
[filled triangle], Cowie et al. (1999) [filled square], 
Gallego et al. (1996) [filled circle]. 

Fig. 9. — The relationship between the gas density 
per comoving volume of the Universe and z com- 
pared to the model predictions (dotted line). The 



gas densities observed from damped Lya systems, 
corrected for missing absorbers, are from Pei et al. 
(1999; filled triangles) and Rao & Turnshek (2000; 
open circles). 

Fig. 10. — The relationship between the 
[Nil] /Ha ratio and the Hubble type for the "qui- 
escent" galaxies from the sample of Kennicutt 
(1992). The large symbols represent averages. 

Fig. 11. — The relationship between the Ha ex- 
tinction as determined from the Balmer decre- 
ment using an underlying Balmer absorption of 
HPE.W. = 2 A (top) and H/SE.W. = 5 A (bot- 
tom) and the Hubble type for the "quiescent" 
galaxies in the sample of Kennicutt (1992). The 
large symbols represent averages. 

Fig. 12. — The relationship between the Ha to 
UV fiux ratio and the axial ratio. The dashed line 
shows the best fit to the data. 
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Table 1 
The Virgo Sample^ 



mag. limit 


No. objects 


Halpha 


UV 


CO 


HI 


H OT K 


11 


3 


3 


3 


3 


3 


3 


12 


8 


8 


6 


8 


8 


8 


13 


25 


22 


19 


19 


25 


25 


14 


45 


29 


27 


23 


45 


42 


15^ 


54 


32 


33 


23 


54 


51 


16^ 


59 


35 


34 


23 


59 


56 



'^including only objects with HI- deficiency <0.3 
''limited to the "ISO" sample 
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Table 2 

The Coma/A1367 supercluster and Cancer sample'' 



mag. limit 


No. objects 


Halpha 


UV 


CO 


HI 


H OT K 


14 


7 


7 




7 


7 


7 


14.5 


20 


20 


5 


17 


20 


20 


15 


63 


50 


13 


42 


63 


63 


15.7 


174 


117 


34 


66 


174 


174 



'^including only objects with Hl-deficiency <0.3 
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Table 3 

The statistics in different morphological bins. 



type 


No. objects 


No. Ha or UV 


Sa 


26 


18 


Sab 


14 


13 


Sb 


29 


19 


Sbc 


36 


19 


Sc 


66 


50 


Scd 


5 


3 


Sd 


3 


2 


Sdm - Sd/Sm 


2 


2 


Sm 


4 


3 


Im - Im/S 


3 


2 


Pec 


26 


21 


S/BCD - dS/BCD 








Sm/BCD 


2 


2 


Im/BCD 


1 





BCD 


2 


1 


S (dS)... 


12 


5 


dIm/dE 


1 





? 


1 


1 
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Table 4 
Adopted IMF parameters 



IMF slope 


IMF cut-off 




Kuvia^Mupf" log[KHa{a, 


M^p)/Kuv{a,M^p)] 


1.5 


80 


1/1.61x10^2 


1/2.01x10"° 


1.903 


2.5 


40 


1/5.41x1040 


1/3.18x1039 


1.231 


2.5 


80 


1/1.16x10"! 


1/3.54x1039 


1.514 


2.5 


120 


1/1.60x10"! 


1/3.66x1039 


1.640 


3.5 


80 


1/5.53x10^^ 


1/I.67xl03« 


0.520 



'^from Chariot & Fall (1993), for Mio^ = 0.1 M©, in units of (M © yr-i)/(ergs-i) 

^from Chariot, private communication, for Miow = 0.1 M©, in units of (M © yic~^)/{ergs~^A.~^) 
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Table 5 

Adopted Ha and UV extinction parameters 



Type 


A{Ha)f.o. 


A{UV)f.o. 


ki{ty) 


Sa-Scd (Pec, S... at D>30 Mpc) 


1.1 


0.60 


0.56 


Sd-Im-BCD 


0.6 


0.20 


0.00 
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